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We develop a three-temperature model to simulate the time dependence of electron and phonon 
temperatures in high-temperature superconductors displaying strong anistropic electron-phonon 
coupling. This model not only takes the tight-binding band structure into account, but also is 
vahd in superconducting state. Based on this model, we calculate the time-resolved spectral func- 
tion via the double-time Green's functions. We find that the dip-hump structure evolves with the 
time delay. More interestingly, new phononic structures are obtained when the phonons are excited 
by a laser field. This signature may serve as a direct evidence for electron-vibration mode coupling. 

PACS numbers: 74.25. Jb, 74.72.-h, 71.38.-k, 79.60.-i 
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The discovery of high-temperature superconductors 
(HTSC) has raised an important issue on the mechanism 
leading to the formation of Cooper pairs, which is still 
under debate. To address this issue, a number of spec- 
troscopy techniques [ll-0| have been used to study the 
role and nature of bosonic modes, to which electrons are 
strongly coupled. In particular, due to technological ad- 
vances and the improved sample quality, angle-resolved 
photoemission spectroscopy (ARPES) [1] has been used 
to probe details of the energy and momentum struc- 
ture of single-particle excitations via the measurement 
of photoemission intensity. However, different interpre- 
tations of the same data may lead to completely different 
mechanisms 0, S]- For example, the dip- hump structure 
observed in ARPES could be interpreted 14| as 

naturnally occurring in the interacting system but having 
little effect on the superconducting pairing mechanism. 
It could also be interpreted [l5| in terms of phonon modes 
that could drive c?-wave pairing. Recent theoretical anal- 
ysis [16!] on the spectral function of thermally excited 
electrons in the cuprates has shown that the out-of-plane 
and out-of-phase buckling mode strongly couples to the 
electronic states near the anti-nodal M points in the Bril- 
louin zone, while the in-plane breathing mode couples 
strongly to the electronic states near the d-wave nodal 
points. The signature of the anistropic electron-phonon 
(el-ph) coupling seems to get enhanced in the supercon- 
ducting state, suggesting the significant role of the el-ph 
interaction in the superconducting mechanism. 

The information extracted from the conventional 
ARPES is limited. Time-resolved ARPES offers the 
capability to simultaneously capture the single-particle 
(frequency domain) and collective (time domain) infor- 
mation, thus making it possible to directly probe the link 
between the collective modes and single-particle states. 
In this setting, either electrons or lattice vibrational 
modes can be selectively excited with an ultrafast laser 
pulse. Recent applications of this technique include the 
studies of transient electronic structure in Mott insula- 
and high-Tc cuprates [l9| with optical pump. 



Furthermore, direct pumping of vibrational mode has 
also been realized in manganites (20| . though not yet in 
the cuprates. Motivated by the thrust of this experi- 
mental technique, in this Letter, we aim to provide a 
theoretical underpinning of transient electronic structure 
for HTSC, in a hope to better understand the nature of 
bosonic modes in these systems. As such, the time evolu- 
tion of the el-ph coupling is investigated for both normal 
and superconducting states with the time-resolved spec- 
tral function. Our theory consists of two parts. First, we 
develop a three-temperature model to simulate the time 
dependence of the electron and phonon temperatures. 
Then, based on the thrce-tempcraturc model, we calcu- 
late the time-resolved spectral function. Our results show 
a kink-structure in the time dependence of the electronic 
temperature at the superconducting transition temper- 
ature. Accordingly, we find from the spectral density 
that the energy position of the phonon mode is offset by 
a time-dependent gap function. More interestingly, new 
signatures of the el-ph coupling, which are absent when 
electrons are excited, can be observed in the case of se- 
lective excitation of phonons. This signature may serve 
as a direct evidence for the electron-vibration mode cou- 
pling as opposed to the coupling between electrons and 
spin fluctuations. 

Three-temperature model: Consider a two-dimensional 
superconductor exposed to a laser field. The model 
Hamiltonian for a vibrational mode v can be written as 
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hilation operators for an electron with momentum k and 
spin a (phonon with momentum q and vibrational mode 
i^), A^q — b'L^q + b^q, the quantity is the normal-state 
energy dispersion, fj, the chemical potential, the gap 



2 



function, and gv the coupling matrix. 

By performing the BogoHubov-de Gennes transforma- 
tion [2l|, Ckt = Ukak - vi^pI. and c_ki = -"kA 
we obtain E, 



Ek £^k((4«k) - (/3k/3,i)), where E^, = 

A^. The time evolution of {a^a-k} and (/3k/?k) 
calculated using the Heisenberg equations-of-motion ap- 
proach. They are found to be 
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and a(^k4)/5T - -9(aj,ak)/9r, where A = /(i;k) = 
+ 1), TVj,^ ^ 7v(r!o) = i/(e^P'^^^o - 1), 5^ = 

6{Eu-q ~ Ek + ilo), and S2 = (5(i?k-q - - ^o)- Here 
we have specifically set il^ ^ Hq. Differentiation of both 
sides of the expression for E^ and substitution of Eq. ^ 
leads to the rate of energy exchange 



dEe 
dr 



in 



no 



5^(ukWk-q - WkWk-q)^'5(£:k-q " -^k - f^o) 



1 



/k(l - /k-q)A^no • 



(3) 



Recent experiment [19] shows that there exists cold lat- 
tice, which is negligibly coupled to the electrons but dis- 
sipates the energy of hot phonons via anharmonic cool- 
ing. Considering this observation, we write a set of rate 
equations for the electron, hot phonon, and lattice tem- 
peratures as 
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where Pe is the power for pumping electrons and Pph the 
power for pumping hot phonons. The specific heat of 
electrons can be calculated from the Boltzmann entropy 

= -2fcBEk{[l-/(^^k)]ln[l-/(£;k)] + /(£;k)ln/(£;k)} 
by Ce — T^dSe/dTf., while the specific heat of hot 
phonons for one- vibrational mode can be calculated from 
the simple relationship Cph — hil.dN{V,)/dTph\n=no- 
The results are given by 
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respectively, with coth x = (e^ + e ^)/ (e^ — e ^). Equa- 
tions constitute our three-temperature model. 
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FIG. 1: (Color) Time evolution of electron (Te), phonon 
{Tph), and lattice (Ti) temperatures for selectively exciting 
electrons (a) and phonons (b) with a laser pulse. 



The original version of this model was phenomeno- 
logically proposed J19j as an extension of the two- 
temperature model [22] for the normal state. The present 
model has the following advantages: (i) it incorporates 
the detailed band structure, (ii) it is valid for both nor- 
mal and superconducting states, and (iii) it includes the 
anistropic effect on the el-ph coupling. 

Now we apply our three-temperature model to sim- 
ulate the time evolution of the electron, hot phonon, 
and lattice temperatures for a d-wave superconductor. 
We use a five-parameter tight-binding model [2^ to 
describe the energy dispersion: ^k = — 2ti(cosfca; -I- 
cosfcj,) — 4t2 cos fca; cos fcj, — 2i3(cos2fcj, + cos2fcj,) — 
4i4(cos 2kx cos fcj,-f cos k^ cos 2ky)—At^ cos 2kx cos 2ky — fj,, 
where h = 1, = -0.2749, = 0.0872, U = 0.0938, 

— —0.0857, and fi — —0.8772. Unless specified exphc- 
itly, the energy is measured in units of ti and the time 
is measured in units of h/ti hereafter. (For ti = 150 
meV, it corresponds to 1740 Kelvin in temperature and 
h/ti to 4.4 femtosecond in time.) The d-wave gap func- 
tion has the form = Ao(Te)(cosfcj, — cosky)/2. The 
temperature-dependent part is given by [24] Ao(Te) = 
Aootanh{(7r/z) v/ar(T,/Te - 1)}, where z - Aoo/(fcBr,) 
and tanhx — 1/cothx. In our calculations, we 
set Aoo = 0.2, the critical temperature Tc = 0.06, 
the specific heat jump at Tc is r = ACe/Ce ^ 
1.43, and a = 2/3. In this work, we focus on 
the buckling phonon mode, for which [3: EH, [2^ 
= 5o{[cos2((7,/2)-Kcos2((7,y/2)]/2}-i/2{^^(k)c/.,(k + 
q)cos((7j,/2) - (/)j^(k)(/)j^(k -I- q)cos(ga;/2)} with 
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txyMxMr] ' , taM = -2tiSm{ka / 2) , and txy,k = 

— 4^2sin(fc^ /2)sin(fc„ / 2) . To be consistent with experi- 
ment [la, |l9|], we set 50 = 0.4 and the mode frequency 
rio = 0.3. The relaxation time Tp = 200. The pump 
power is a Gaussian pulse of P = Poe~'^^/'-^°'^\ for which, 
the corresponding FWHM (full width at half maximum) 
is 2.35 a. To pump electrons, we set Pq = 0.15 and 
a — 1. Considering that the energy scale of phonons is 
much smaller than that of electrons, we set Pq — 0.006 
and a = 100 for pumping phonons. The ratio Cph/Ci in 
Eq. dH is set to be 0.2. 
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Figure [T] displays the temporal evolution of three re- 
spective temperatures for pumping electrons (a) and hot 
phonons (b). Before pumping, all three types of de- 
grees of freedom (DoF) are in equilibrium, which is set at 
Te = Tph ^Ti ^ 0.01. As shown in Fig. [^a), when the 
pump pulse is absorbed by the electrons, the electronic 
temperature rises steeply around r = 0, and reaches the 
maximum after a small time delay. It then begins to drop 
at a time scale determined by the el-ph coupling strength, 
followed by a slower relaxation. Interestngly, we also ob- 
serve a small kink in Te at Tg = Tc, which entirely arises 
from the breakup of the Cooper pairs, resulting in the 
dramatic change in the rise rate of Te- This kink was not 
captured in the early model [10] . Simultaneously, the hot 
phononic temperature T^^ rises smoothly via the energy 
exchange with electrons and then drops by exchanging 
energy with the cold lattice, causing the slight increase 
of Ti. When directly pumping phonons, the time depen- 
dence of the temperatures is similarly observed, including 
the kink in T^, as shown in Fig. [TJb). Due to the large 
pumping width, the kink is more easily visible in this 
case. 

Time-resolved spectral function: The time-resolved 
spectral function is defined as A(k, = — ^ Im^n (k, w) . 
Here Qn is the one-one component of the retarded 
Green's function ^(k, cj), which is related to the self- 
energy by ^"^(k, cj) = ^(7\k, - E(k, cj), with 

^(^"^(k,^) = wCTo - AkiTi - CkCTs- o'o, 1,2,3 are the unit 
and Pauli matrices. As is known [l6| . in the equilibrium 
state with Tg = Tp^, the self-energy can be evaluated 
more conveniently within the imaginary-time Green's 
function approach, in which a key step is to convert 
the Bose-Einstein distribution to the Fermi distribution, 
nB{iiOn±Eif,^q) = -"'F(±£^k-q) (with w„ = {2n + l)'KT, 
T = = Tph). For the current situation, the temper- 
atures of electrons and hot phonons are no longer tied 
to each other and the above conversions are not valid 
any more. To avoid this restriction, here we apply the 
double-time Green's function approach ^\ to calculate 
S(k,w). The result is 

^ q 

+[2i;k_q($i - $3) + 2no(*3 - *4)]cti 

+[ek($i - $2) + (f^oek/i^k)($3 - $4)1^3}, (9) 

where $1 = 7V(f]o)/[(a; - Oq + E^^^){lo - Oq - £^k-q)], 
$2 = iV(-r!o)/[(^ + f^o + £^k-q)(w + r!o - i^k-q)], *3 = 

/(-£;k_q)/[(cj + - £^k-q)(w - l^O - Sk-q)] , aud $4 = 
/(i;k-q)/[(w + + £^k~q)(w - Oo + i^k-q)]. 

For simplicity, let us first look at the density of states 
(DOS), which is defined by p{uj) = ^ Y.\^ A{\<i, oj) . For 
direct pumping of electrons, we calculate the DOS at time 
sequences r = —100, —2.7, 2.8, and 600, respectively. 
The results are plotted in Fig. [5] (a) . At the initial time 
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FIG. 2: (Color) Time evolution of the density of states for se- 
lectively exciting electrons (a) and phonons (b), respectively. 



T — —50 where = 0.01 and the system is at supercon- 
ducting state [Ao(Te) — 0.19], we observe two signatures 
of the el-ph coupling located at ±(rio+Ao(Te)), as shown 
in Fig. [2] (dip 1 and dip 2). They are symmetric with re- 
spect to cl) = 0, but dip 2 is much stronger. It is related 
to the fact that the normal-state van Hove singularity is 
located below the Fermi energy and as such the coher- 
ent peak at — Ao(Te) has stronger intensity than that at 
Ao(Te). In addition, dip 3 occurs near the characteris- 
tic energy —{Em + f^o)j where Em is the quasiparticle 
energy at k = (7r,0) or equivalent wavevector points in 
Brillouin zone. It arises from the van Hove singularity at 
k = (tt, 0). The signature at the energy Em + f^o is much 
weaker also because of the van Hove singularity location 
in the normal state. At r = —2.7 where Tp = 0.041 and 
Ao(T'e) ~ 0.1, similar behaviors are observed. However, 
at r = 2.8, Tg = 1.05 and thus only the signature of the 
normal-state el-ph coupling (dip 3) can be observed. At 
T = 600, a similar structure of the normal-state el-ph 
coupling is observed. 

For direct pumping of phonons, we choose different 
time sequences t = —1000, —158, 110, and 800 to eval- 
uate the DOS. The results are displayed in Fig. [H^b). 
From Fig. [1Kb), we observe that the dips suggesting 
the el-ph couplig in both normal and supercondcting 
states becomes stronger than those for selectively pump- 
ing electrons, as expected. In particular, we observe two 
new small dips aX lu — ±0.05 for r = —158 (at which 
Te = 0.041 and Tph = 0.163), both of which arise from 
the poles of Sk(w) at w = ±(r2o — Ao(re)), in addition to 
those observed in the case of electrons being excited di- 
rectly. This is because the hot phononic temperature be- 
comes very high while the electronic temperature is still 
cold and the contributions from the terms as weighted by 
A^(±ilo) are significantly enhanced. Figure. H^b) clearly 
shows that, as time elapses from r = —1000, —158, to 
110, the peaks move toward the zero energy, while from 
T = 110 to 800, the locations of these two peaks remain 
unchanged. This can be understood by considering that, 
after pumping, T^ rises with the time delay, while Ao(Te) 
decreases and rapidly vanishes when Tg reaches Tc- 

Finally we numerically evaluate the time-resolved spec- 
tral function for the two excitations discussed above at 
those time sequences chosen for calculating the DOS. In 
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FIG. 3: (Color) Time-resolved spectral function for selectively 
pumping electrons (ai-a4 for t = —100,-2.7,2.8,600) and 
phonons (bi-b4 for r — —1000,-158,110,800), respectively. 
The vertical axis is kx/ir at a fixed ky = O.TStt and the hori- 
zontal axis is energy. The black, red, and magenta lines rep- 
resent the energy location, —{Em + ^1o), — (Ao(Te) + r2o), and 
— {flo — Ao(re)), corresponding to the dip or hump location 
discussed in Fig. [2] 



our calculations, we consider cuts along A;a;-axis in the 
Brillouin zone at a ky chosen near the zone boundary. 
Fig. [3] shows the snapshots of the image of the spec- 
tral function A(k,a;) as a function of and u. From 
Fig. [31[ai,a2), we observe kinks at uj — — (Ao(Te) 4- r^o) 
(red lines) and —{Em + ^o) (black lines) at r = —100 
and —2.7. The red line is shifted in energy as time elapses 
from r = -100 to r = -2.7 because the BCS gap Ao(re) 
drops significantly from 0.19 to 0.1. These kinks cor- 
respond to the dip structures discussed in Fig. [51 As 
T > —2.3, the superconducting gap vanishes, so do the 
kinks marked with red lines (see Fig.[3l^a3-a4)). The kinks 
marked with black lines persist through the whole time 
sequence, which can be ascribed to the el-ph coupling 
signature in the normal state. Therefore we have ev- 
ery reason to regard these kinks marked with red lines 
as the signature of the el-ph coupling in the supercon- 



ducting state. For hot phonon pumping (Fig. [31^bi-b4)), 
in addition to those signatures of the el-ph coupling 
observed for electron pumping, new kink structure at 
— (rjg — Ao(Te)) (marked with magenta line in Fig.[2)Jb2) 
is also observed. These observations are consistent with 
our finding in the DOS. In the viewpoint that the elec- 
tronic spin fluctuations arise from the strong correlation 
within the electronic DoF itself, and if one can assume 
that the spin fluctuations will ride on the electrons and 
thus its effective temperature will be tied to the electronic 
temperature, direct pumping of hot phonons will provide 
a unique way to differentiating the bosonic modes be- 
ing of electronic or phononic origin, to which electronic 
quasiparticles are strongly coupled in HTSC. 

In conclusion, we have developed a three-temperature 
model to simulate the time evolution of temperature for 
electrons, hot phonons, and the lattice. Our model takes 
both quasiparticle excitation and relaxation into account 
and thus is valid in both normal and superconducting 
states, because at the early stage of excitation, the sys- 
tem can be still in the supercondcting state. Based on 
this model, we derive a formula for the time-resolved 
spectral function, allowing us to investigate the dynamics 
of the el-ph coupling for HTSC materials. 
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